β-BaB2O4 with special cut-angle applied to single crystal cascaded third-harmonic generation
Ren Hong-Kai1, 2, Qi Hong-Wei3, Wang Zheng-Ping1, 2, †, Wu Zhi-Xin1, 2, Wang Meng-Xia1, 2, Sun Yu-Xiang1, 2, Sun Xun1, 2, Xu Xin-Guang1, 2, ‡
State Key Laboratory of Crystal Materials, Shandong University, Jinan 250100, China
Key Laboratory of Functional Crystal Materials and Device (Shandong University), Ministry of Education, Jinan 250100, China
School of Chemistry and Chemical Engineering, Linyi University, Linyi 276000, China

 

† Corresponding author. E-mail: zpwang@sdu.edu.cn xgxu@sdu.edu.cn

Project supported by the Natural Science Foundation of Shandong Province, China (Grant Nos. ZR2017MF031 and ZR2018BF029).

Abstract

High-efficiency single crystal cascaded third-harmonic generation (THG) was realized in β-BaB2O4 (BBO) material with special cut-angle. By analyzing effective nonlinear optical coefficient (deff) of the cascaded THG process, which was composed by type-II frequency doubling and type-I sum-frequency, the optimum phase matching (PM) direction in BBO crystal was determined to be (θ = 32.1°, ϕ = 11°). With an optimized 9-mm long sample which was processed along this direction, the highest cascaded THG conversion efficiency reached 42.3%, which is much superior to the similar components reported previously, including ADP, KDP, and GdxY1−xCOB crystals.

1. Introduction

Now, the bulk nonlinear crystals are still the primary mediums for optical frequency up-conversions. The popular nonlinear optical (NLO) crystals include KDP, KTP, β-BaB2O4 (BBO), LBO, and so on. For the third-harmonic generation (THG) of near infrared laser, traditionally two NLO crystals are necessary: the first is for second-harmonic generation (SHG) and the second is for sum-frequency generation (SFG).[16] Recently, we developed a cascaded THG apparatus that only uses one NLO crystal. It possesses some advantages such as wide available waveband, high efficiency, and cost saving. At the same time, it has been successfully applied to several NLO materials, including KDP, ADP, and GdxY1−xCOB.[79] In present paper, for this design the BBO crystal was used as the NLO medium which brought much larger nonlinearity, and at the same time a special coated quarter-wave plate (QWP) was employedto permit completely utilizing the fundamental wave. As a result, with the optimized BBO sample, we achieved a highest conversion efficiency of 42.3%, which is much better than those previously obtained from other NLO materials.

2. Theoretical analysis of BBO crystal
2.1. Phase-matching angle

According to the crystallography tradition, the polar coordinates (θ, ϕ) are used to denote the crystal spatial direction, where θ is the included angle with the Z-axis, and ϕ is the azimuthal angle in the plane perpendicular to the Z-axis. As a trigonal symmetric, optical uniaxial crystal, the Z-axis of BBO crystal is coincident with the optical axis, and the phase matching (PM) condition is only determined by the parameter θ. Based on the Sellmeier equations of BBO crystal,[10] the PM angles θ are calculated as a function of fundament wavelengths for the SHG (ω + ω → 2ω) and SFG (ω + 2ω → 3ω) processes respectively, as shown in Fig. 1. Among all of the five PM styles, type-II SHG and type-I THG present the closest PM directions, as indicated by the red and blue lines in Fig. 1. The discrepancy of their PM angles θ is −3.0° at 750 nm and 2.9° at 1200 nm, respectively. If Δθ < 3° is used as an application standard of one crystal cascaded THG technique, the work bandwidth of BBO crystal is larger than 450 nm, which covers 750–1200 nm. This value is much better than that of the KDP or ADP crystal, whose work bandwidth is about 100 nm, i.e., the available fundamental wavelength is from 1000 to 1100 nm.[7,8] At the wavelength of Nd:YAG laser (1064 nm), the θ value of BBO crystal is 32.9° for type-II SHG and 31.3° for type-I SFG. It means that the θ direction of the practical processed sample should be controlled between these two angles.

Fig. 1. (color online) Variation of phase-matching angle with the fundamental wavelength of BBO crystal.
2.2. Effective nonlinear optical coefficient

For type-II SHG and type-I SFG, the effective nonlinear optical coefficients (deff) of BBO crystal can be calculated from[11] where d22 = 2.3 pm/V and d31 = −0.16 pm/V.[12,13] When the θ angle is fixed, the largest deff(II) is at ϕ = 0° direction, while the largest deff(I) is at ϕ = 30° direction. Therefore, an optimized ϕ angle should be chosen to realize the largest THG output for the present one BBO crystal, cascaded frequency upconversions. Such a problem has never been met in previous investigations of KDP, ADP, and GxY1−xCOB crystals, because for those crystals the largest deff values for SHG and SFG appear at the same spatial PM angle.[79] In a plane-wave fixed-field approximation, the SHG output E2 and SFG output E3 have the disciplines where E1 is the original incident fundamental energy, and is the residual fundamental energy after the SHG process, which attends the followed SFG process. Under the ideal THG experimental conditions, .[14] Taking this relationship and Eq. (3) into Eq. (4), there is Utilizing Eqs. (1) and (2), the variation of (deff, I)2(deff, II)4 with ϕ angle (0–30°) can be obtained, and the result is shown in Fig. 2. In this calculation, the θ angle is fixed to be 32.1°. It can be seen that for the highest THG output the largest (deff, I)2(deff, II)4 value appears at ϕ = 11°, not the middle angle ϕ = 15°. At (θ = 32.1°, ϕ = 11°) direction, the |deff, II| and |deff, I| values are 1.38 and 0.79 pm/V respectively, which are much larger than those of KDP, ADP, and GdxY1−xCOB crystals (< 0.6 pm/V).

Fig. 2. (color online) Variation of (deff, I)2(deff, II)4 with the ϕ angle of BBO crystal.
2.3. Effective lengths[15]

The quasi-static interaction length Lqs can be calculated from where τ is the radiation pulse duration and ν is the inverse group-velocity mismatch. For BBO crystal, the parameter ν is in the order of 102–103 fs/mm, so under our experimental condition, i.e., 40 ps of FWHM (full width at half maximum) pulse duration, Lqs is 4–40 cm. The actual sample length used in our experiments is no more than 1.5 cm, so the present frequency conversion processes coincide with the quasi-static interaction condition, and the group-velocity mismatch can be ignored.

The aperture length La can be obtained from where d0 is the beam diameter and ρ is the “walk-off” angle. For type-II SHG and subsequent type-I SFG of 1 μm laser, ρ is about 4° for BBO crystal. In our experiments, d0 is 2 mm, so La is calculated to be 28.6 mm and the practical sample should be shorter than this length.

For the nonlinear interaction length LNL, there is where σ3 is the nonlinear coupling coefficient of the generated wave which is proportional to the effective nonlinear optical coefficients deff, and inversely proportional to the refractive index n. For KDP crystal (deff ≈ 0.35 pm/V, n ≈ 1.5), σ3 is about 10−6V−1. Therefore, for the (θ = 32.1°, ϕ = 11°) oriented BBO crystal discussed above (deff ≈ 0.8–1.4 pm/V, n ≈ 1.6), σ3 will be 2 × 10−6–4 × 10−6 V−1. In Eq. (8), α0 is the wave amplitude of fundamental laser at the input surface of the crystal, which can be calculated by where P1 is the power of the fundamental radiation and ω is the beam radius. Taking Eq. (9) into Eq. (8), the variation of LNL as a function of the input power density P1/ can be plotted, as shown in Fig. 3. In this calculation about BBO crystal, we choose an average σ3 value of 3 × 10−6 V−1 and the n value of 1.6. In our experiments, the highest crystal incident energy is 4 mJ at a full temporal duration of 100 ps and a beam radius of 1 mm, so the largest input power density is 1.3 × 109 W/cm2. From Fig. 3, it can be seen that the corresponding LNL is 4.3 mm, which means that in theory the optimum crystal length will be 8.6 mm or so (2LNL).

Fig. 3. (color online) Variation of the nonlinear interaction length LNL with the fundamental power density of BBO crystal.
3. Experiment and discussion

The experimental equipment is similar with that used for KDP crystal. Besides different NLO crystal, an important change occurs on the 532 nm quarter wave plate (QWP). Previously, both the fundamental wave (1064 nm) and the 2ω wave (532 nm) passed through the QWP to and fro, which led to the fact that the polarization direction of fundamental wave cannot be optimized completely. Now its surface towards the NLO crystal is high-reflective (HR) coated at 1064 nm and anti-reflective (AR) coated at 532 nm, and the other surface towards the total reflection mirror M2 is AR coated at 532 nm. The fundamental wave no longer passes though the 2ω QWP. In this way, the incident fundamental polarization can be adjusted freely, and the unfavorable fundamental polarization variation induced by the 2ω QWP is avoided, i.e., the optimal scheme (OS) is realized.[7] At the same time, the 2ω QWP is tightly attached to the mirror M2 to decrease the optical path difference between the ω and 2ω waves.

We prepared six BBO samples for this experiment, which were cut along two directions, (32.1°, 11°) and (32.1°, 15°). For each direction, three crystal lengths were selected, 6 mm, 9 mm, and 12 mm. The cascaded THG results are shown in Fig. 4. All of the six samples exhibit efficient cascaded THG conversion with conversion efficiencies larger than 10%, which proves that the selection of θ angle from the Fig. 1 is reliable. Generally, at the same crystal length, the (32.1°, 11°) sample always presents better performance than the (32.1°, 15°) sample. For example, the largest THG output and the highest conversion efficiency of the 6 mm, (32.1°, 11°) BBO are 0.95 mJ and 38.1%, respectively, and the corresponding data of the 6 mm, (32.1°, 15°) BBO are 0.62 mJ and 20.0%. Such property proves the theoretical analysis result of Fig. 2, i.e., the (32.1°, 11°) direction possesses larger combined NLO coefficient than the (32.1°, 15°) direction for the cascaded THG process. For BBO samples at the same direction, the output energy and conversion efficiency of the 9 mm crystal are obviously superior to the results of 6 mm and 12 mm crystals. It experimentally manifests that the optimized crystal length is 9 mm or so, which is in good agreement with Fig. 3. To sum up, the best experiment data come from the (32.1°, 11°), 9 mm sample, i.e., the largest THG output of 1.17 mJ at a fundamental energy of 3.51 mJ, and the highest conversion efficiency of 42.3% at a fundamental energy of 1.02 mJ. The most inferior are from the (32.1°, 15°), 12 mm sample.

Fig. 4. (color online) Variation of (a) cascaded THG output energy and (b) conversion efficiency with the fundamental energy of different BBO samples.
4. Conclusion

BBO crystal is one of the most popular ultra-violet NLO crystals. In this paper, we realized its single block, cascaded THG for the first time. By optimizing the crystal PM direction and the crystal length, as well as utilizing the polarization OS, the large nonlinearity of BBO crystal is fully exploited and a high THG conversion efficiency of 42.3% is obtained. Comparing with the traditional cascaded THG converter composed by two NLO crystals, the present design exhibits some advantages such as compact, cost saving, and super wide applicable waveband (>450 nm). Its work wavelength extends to 800 nm around, so efficient single crystal, THG converter maybe developed for ultrafast Ti:sapphire or diode lasers in the future. Besides the practical value, this study will also lend good design reference for other cascaded frequency upconversion processes and NLO materials.

Reference
[1] Yu Z H Tian J R Song Y R 2014 Chin. Phys. 23 094206
[2] Liu G Q Ying H P 2014 Chin. Phys. 23 050502
[3] Zhao S M Zhuang P 2014 Chin. Phys. 23 054203
[4] Ueda K Orii Y Takahashi Y 2016 Opt. Express 24 30465
[5] Feng X Zhu L L Li F Q 2016 RSC Adv. 6 33983
[6] Chen J Zheng Y L An N 2015 Opt. Lett. 40 4484
[7] Zhang L Li K Xu D G 2014 Opt. Express 22 17187
[8] Gapontsev V P Tyrtyshnyy V A Vershinin O I 2013 Opt. Express 21 3715
[9] Furuya H Nakao H Yamada I Ruan Y F 2000 Opt. Lett. 25 1588
[10] Qi H W Wang Z P Yu F P 2016 Opt. Lett. 41 5823
[11] Qi H W Sun Y X Wang Z P 2017 Chin. Phys. 26 084201
[12] Qi H W Wang Z P Yu F P 2016 IEEE Photonics J. 8 1602107
[13] Kato K 1986 IEEE J. Quantum Electron. QE-22 1013
[14] Midwinter J E Warner J 1965 Brit. J.Appl. Phys. 16 1135
[15] Roberts D A 1992 IEEE J. Quantum Electron. 28 2057
[16] Chen C T Wu B Jiang A 1985 Sci. Sinica Ser. 28 235
[17] Craxton R S 1981 IEEE J. Quantum Electron. 17 1771
[18] Dmitriev V G Gurzadyan G G Nikogosyan D N 1999 Handbook of Nonlinear Optical Crystals III Berlin Springer 10.1109/JQE.1981.1071318